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Abstract 

We analyse parafcrmionic operators in the Q -state Potts model from three different perspec- 
tives. First, we explicitly construct lattice holomorphic observables in the Fortuin-Kasteleyn 
representation, and point out some special simplifying features of the particular case Q = 2 
(Ising model). In particular, away from criticality, we find a lattice generalisation of the 
massive Majorana fermion equation. We also compare the parafcrmionic scaling dimensions 
with known results from CFT and Coulomb gas methods in the continuum. Finally, we 
show that expectation values of these parafermions correspond to local observables of the 
SLE process which is conjectured to describe the scaling limit of the Q-state Potts model. 
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1 Introduction 



Among the special features of two-dimensional physics, one of the most intriguing is semi- 
locality, which is marked by the appearance of a complex phase in the correlation function of 
two or more excitations once their positions have been exchanged. The study of this type of 
excitations, called parafermions and characterized by fractional spin, has found fruitful applica- 
tions in diverse areas of theoretical physics. Parafermions are rather well understood in critical 
systems in the continuum limit, since the discrete symmetries they are associated to combine 
with conformal invariance and permit an extensive classification of the so-called parafermionic 
Conformal Field Theories. The first of such theories, constructed by Fateev and Zamolodchikov 
m, describes the critical point of certain lattice models with Z^r symmetry, which reduce to 
the Ising and 3-state Potts models for N = 2 and = 3, respectively. The parafermionic 
fields represent the chiral conserved currents associated to the extra symmetry, the conservation 
laws being expressed in complex coordinates as a holomorphicity condition on the parafermions. 
In the Ising case, the parafermionic current is precisely the chiral component of the Majorana 
fermion that is known to describe the model at criticality and in the neighbouring scaling regime, 
and the holomorphicity condition is deformed, away from criticality, by terms proportional to 
the fermion mass. 

In this paper, we analyse parafermionic observables in the Q-state Potts model, both on the 
lattice and in the continuum limit. With respect to the literature on the subject, our approach 
provides a new explicit proof of holomorphicity on the lattice, and the first analysis of its relation 
to Stochastic Loewner Evolution (SLE) in the continuum limit. 

On the lattice, we adopt the Fortuin-Kasteleyen (FK) representation and define the parafermions 
in terms of the corresponding gas of loops. This allows us to show that, for every value of Q, 
there is a particular value of the spin at which the parafermionic correlation functions are purely 
holomorphic (in a lattice sense) at criticality. The resulting scaling dimensions agree with CFT 
and Coulomb gas results in the continuum. 

It is conjectured that the continuum limit of the Q-state Potts model in the FK representa- 
tion is described by SLE^ with ^/Q = —2 cos (Att/k), 4 < k < 8. We provide further evidence 
for this correspondence by constructing SLE observables which are holomorphic, have the appro- 
priate scaling dimension and satisfy the same boundary conditions as the lattice parafermions 
mentioned above. Finally, we explicitly discuss the relation between parafermions and SLE on 
the lattice in a particular example. 

Our work, although carried out independently, has some overlap with that of the programme 
of Smirnov[2] for proving that the scaling limit of certain curves in lattice 0(n) and Q-state Potts 
models is SLE with the appropriate value of k. (which has been pushed to completion by Smirnov 
for the Ising case {n = 1, Q = 2)). Our emphasis, however, is on the lattice holomorphicity of the 
multi-point correlations of the parafermions, and their relations to earlier physics CFT literature, 
rather than on detailed considerations of the one-point function in an arbitrary domain which 
are required to prove convergence to SLE. That fact that analogous holomorphic observables 
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exist in SLE was already pointed out by B. Doyon and the present authors in an earher paper [3]. 

The paper is organized as follows. In Section[2] we briefly review the FK representation of 
the Potts model. In Section[3]we define the parafermionic operators of spin p, and we show that 
their correlation functions are purely holomorphic at a particular value of p. Special features 
of the Ising model are discussed in Section^ Section[5] presents the continuum results, flrst the 
ones already known from Coulomb gas and CFT methods, then the new prediction obtained by 
SLE. Finally, in Section[6] we discuss the relation between the SLE and lattice results, and we 
conclude the paper in Section[71 

2 Q-state Potts model 

Let us consider the Q-state Potts model on the square lattice, at each site of which is associated 
a spin variable Si S {0, 1, Q — 1}. The partition function is 

Z = Y^e<^'^> = JJ + , u = e^-l. (2.1) 

{5} {5} <id> 

By expanding the product in (j2.ip . one obtains the Fortuin-Kasteleyn (FK) representation [1] 

Z = ^u^Q^ (2.2) 

G 

where G is any subgraph of the original domain, consisting of all the sites and some bonds placed 
arbitrarily on the lattice edges, h is the number of bonds in G, and c is the number of clusters 
of connected sites into which the bonds partition the lattice (FK clusters, which also include 
single sites). Fig. [1] shows a configuration in this expansion, where wired boundary conditions 
have been assigned (i.e. n = co on the boundary edges). Alternatively, one could assign n = 
to the boundary edges, which corresponds to free boundary conditions on the Potts spins. 

Although the original formulation of the model requires Q to be a positive integer, the FK 
representation (|2.2p allows one to interpret Q as taking arbitrary complex values. The model is 
known to have a critical point with a diverging correlation length at the self-dual value of u for 
Q real and < Q < 4. 

The partition function (j2.2|) can be also expressed as that of a gas of fully packed loops on 
the medial lattice (the dashed lines in Fig.[T]), whose number is d = / + c for each given G, where 
/ is the number of independent loops in G. 

By using Euler's relation 

c-l = N -b , (2.3) 

where N is the number of vertices in G (i.e. the total number of vertices in the lattice), the 
partition function can be rewritten as 
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Figure 1: A configuration in the FK representation of the Potts model with wired boundary 
conditions (thick hues). Dots and crosses represent sites of the original and dual lattice, respec- 
tively. 



By reexpressing ()2.4p in terms of the dual lattice and imposing self duality at the critical point, 
one extracts the critical value 



which implies that the weights for the FK clusters at criticality are equivalent to counting each 
loop on the medial lattice with a fugacity \fQ. 

3 Holomorphic observables on the lattice 

In this Section, we will recall the definition of parafermion (or spinor) operators, and we will give 
it an explicit geometrical meaning in the context of the FK representation. We will then show 
that, at a particular value of the spin, the parafermions are purely holomorphic observables in 
a lattice sense. 

Let us first briefly recall the definitions of order and disorder operators, which are the basic 
ingredients in the construction of lattice parafermions (for more details, see [SJE])- An order 
(or spin) operator Oi is associated to every site of the original lattice (dots in Fig.[T]), taking 
e Q L where Si e {0, - 1}. The spin-spin correlation function (cTjO"*) = 

{e ^ ) is proportional to a sum like ()2.2p where only those graphs G C {G} appear 

on which the sites i and j belong to the same FK cluster. The contributions of all other 
terms vanish, since Si and Sj belong to different clusters and are freely summed over. Disorder 
operators /i^ can be associated to the sites of the dual lattice (crosses in Fig.[T]) by defining 
their correlation function as (/i^^^*) = Z' /Z , where Z' is obtained from (j2.1|) by modifying 
6{Sk,Si) — > 5{Sk + 1 {mod Q), Si) along a path on the dual lattice that connects sites k and 
/. It has been shown [5l [6] that this definition of the correlation function does not depend on 
the particular path chosen, and that the resulting operator n is dual to the spin operator a. 
Hence the disorder correlation function only gets contributions from those graphs where the 




(2.5) 



sites belong to the same cluster on the dual lattice. 

If we now look at the covering lattice, which consists of the sites of both the original and the 
dual lattice, we can define parafermion (or spinor) operators ipp{e) on each of its edges e as the 
product of the spin and the disorder operator at the sites connected by e, times a phase factor 
which encodes p: 

ijpie) = a{e)fi{e)e-'P'^'''''K (3.1) 

The meaning of the angle ^(7, e) will be clarified below. The existence of several spinor operators, 
labelled by p, at each value of Q, is due to the fact that the disorder operator can be modified by 
replacing 6{Sk, Si) with 6{Sk + n (modQ), Si), with n = 1, 2, Q. It was shown in [6] that the 
different parafermions can be labelled by a fractional index p which has the physical meaning of 
spin, since the spinor correlation function acquires a phase factor e*^'^^ when the spinors circle 
one another. 

In order to give a precise meaning to the angle ^(7, e) in (j3.ip . we now geometrically identify 
the spinor operators in terms of the fully packed loops in the FK representation of the Potts 
model (dashed lines in Fig.[T]). It follows from (jS.ip that spinor correlation functions are observ- 
ables of such loops, since they only get contributions from terms in ()2.2p where all order and all 
disorder operators belong to the same FK cluster on the original and dual lattice, respectively, 
hence all edges e belong to the same loop on the medial lattice. For each term in i\2.2\\ which 
contributes to a spinor correlation function, we can define the parafermions' phases as 

arg[V'p(e)] = -p9{^,e) , (3.2) 

where 7 indicates the loop which crosses all edges involved in the given correlation function, and 
the angle ^(7, e) is the inclination at which 7 intersects the edge e. In the cases of interest p 
is fractional and (13. 2p defines the left hand side only up to a multiple of 27rp. However we can 
give it a unique meaning as follows: for a curve 7 which begins and ends on the boundary we 
choose 0(7, e) for e on the boundary to be in the range [0, 27r), and then define it recursively at 
other points by requiring that between neighbouring edges on the medial lattice it can change 
only by ib7r/2. With this definition we see that 9{'j,e) can be arbitrarily large, depending on 
the degree of winding of 7 around the edge e. If 7 is a closed loop, we choose a given edge e at 
random, require again that ^(7, e) € [0, 2tt), and define the angles of the other edges recursively 
as before. Since all observables will turn out to depend only on the differences of the 9{'y,e) 
between different edges on the same curve, this arbitrariness is not important. 

3.1 Lattice Holomorphicity 

We now show that, at the critical point, ('i/'p(ei)'0p(c2) ■ ■ ■ V'p(cn)) is a "lattice-holomorphic" 
function of each of its n variables, i.e. it satisfies a discrete version of the Cauchy-Riemann 
equations, if we properly choose the value of p. To this purpose, we will prove that, if 
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then 

^ (V'p(ei) • • • V'p(en) V'p(e)) dz^ = , (3.4) 

eeC 

for every closed contour C of the covering lattice which does not include ei, eJl]. This equation 
has a direct meaning as lattice version of a vanishing contour integral. Furthermore, if we 
consider C as the boundary of an elementary plaquette of the covering lattice, (j3.4p is equivalent 
to the discrete Cauchy-Riemann equations, which constitute its real and imaginary part (for a 
systematic study of discrete Cauchy-Riemann equations see [7]). 

Proof of l^.^p . Let us first consider the case when C is the boundary of an elementary 
plaquette of the covering lattice. A loop 7, which intersects all edges ei, ...,e„ and e, can visit 
the plaquette in eight inequivalent ways, which are depicted in Fig. [21 where 7 is represented 
by the continuous lin^. The dashed lines represent other loops in the same configuration; their 
paths outside the plaquette can be different from the picture, but their particular shape does 
not matter, as we shall see below. To each edge it is associated the value of the inclination at 
which 7 intersects the edge in the considered configuration, up to a multiple of 27r, which is, 
however, the same for every edge in the plaquette which the curve visits. 




Figure 2: Inequivalent configurations on a single plaquette 

We will now show that configurations (la) and (lb) cancel each other in the sum l\3A\\ . and 
the same can be verified in the other cases. The result automatically holds in all winding classes, 
since winding simply implies a global rotation of all angles. 

^Singularities occur in the correlation function when two of its arguments coincide. 

^For definiteness, we have oriented the loop so the sites of the original lattice, carrying the Potts spins, lie 
immediately to its left, and the dual sites to its immediate right. 
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The sum over (la) and (lb) gives 



X] (V'p(ei) • • • ^p(en) ^p(e)) fee 



(3.5) 



laAb 



(l + ie'P^'^ P{la) + (l + «e*P5 + (-l)e^P^ + (-i)e-^P5) P(16) 



where -P(la) and ^^(16) are the weights of all graphs in the class (la) or (16), respectively. From 
()2.4|2.5p we know that, at criticality, 



P{la) = y'QPilb) , (3.6) 

since the two types of configuration differ from the presence/absence of a closed loop, weighted 
with a factor ^/Q (independently of its particular shape). Eq. ()3.6p can be interpreted as an 
"equation of motion", since it expresses the variation of the partition function under a local 
rearrangement of configuration. Therefore, (j3.5p vanishes for y/Q = 2 sin(p^), i.e. when (j3.3p 
holds. It is easy to check that the same cancellation mechanism takes place for the other pairs 
of configurations, which proves (13.40 for the single plaquette case. 

The validity of (j3.4p for an arbitrary closed lattice contour C follows directly since it is equal 
to the sum of the contour integrals around each elementary plaquette contained inside C. 

□ 



3.2 One— point function 

By choosing suitable boundary conditions, it is possible to obtain a non-vanishing expectation 
value of ipp. We will now focus on this situation, in view of its simplicity and its relation to SLE, 
which will be discussed in Section[6l If we impose wired boundary conditions on one connected 
component of the boundary and free boundary conditions on its complement (see Fig. [3]), a 
domain wall 7 will be generated on the medial lattice (continuous line), which is simultaneously 
the boundary of the cluster attached to the wired part of the boundary, and of the dual cluster 
attached to the rest of the boundary. In addition, of course, there will in general be closed loops 
(dashed lines). 

The expectation value {ipp{e)) vanishes every time the values of cr{e) and/or fj.{e) can be 
summed over freely in the partition function. As a consequence, (V'p(e)) is different from zero 
only on the edges crossed by the domain wall, since in this case both a{e) and fi{e) belong to a 
cluster connected to the boundary and therefore have a fixed value (wired boundary conditions 
are equivalent to fixing the value of a and free boundary conditions to fixing fi). Therefore, we 
can define the phase ^(7, e) in (|3.2p as the inclination at which the domain wall 7 intersects the 
edge e. It is easy to check that the proof of (j3.4p also holds for 

(M^)) Sze = 0, (3.7) 
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Figure 3: A configuration in the FK representation of the Potts model with wired and free 
boundary conditions (thick lines on the original and dual lattice, respectively) 

by simply considering 7 as the domain wall described above. Therefore, (^p(e)) is a lattice 
holomorphic function. 

Moreover, we can determine which are the boundary conditions satisfied by {'ipp{e)). To 
this purpose, we have to look at the plaquettes of the covering lattice which are cut in half by 
the boundary edges of the original or dual lattice, and therefore have a triangular shape. The 
presence of the boundary implies that {ipp{e)) only differs by a phase on the two edges of any of 
these plaquettes, since the domain wall automatically intersects the second edge of a plaquette 
once it entered the first. Therefore, we can limit our analysis to half of these edges, and we 
choose for convenience the entering edges of the plaquettes which lie on the left side of the 
boundary (with wired boundary conditions), and the exiting edges cr of the plaquettes which 
lie on the right side (with free boundary conditions). Since winding is not possible around these 
plaquettes, the phase ^(7, e) is uniquely determined by the boundary edge under consideration: 

e{eL) = et-7T, e{eR) = e,, , 

where 9i, is the inclination of the boundary itself (measured counterclockwise from the starting 
point of the curve) and we set to zero the value of 6{'y,e) at the starting edge in Fig.[3] (these 
arguments can be easily generalized to more general domains than the one in Fig. [3]). Therefore, 
the phase of {'ipp{e)) on the boundary is fixed as 

arg[(V'p(eL))] = -p{9b - vr) , arg[(V'p(eR))] = -p9b . (3.8) 

This is the generalization to an arbitrary domain of the boundary conditions satisfied by the 
function z~P on the upper half plane when the domain wall starts from x = on the real axis, 
since in that geometry our convention implies 06 = for x > and 0^ = 27r for x < 0. 
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4 The Ising model 



In this Section, we will analyse in further detail the properties of {ipp{e)), enlightening some 
special features of the case Q = 2, i.e. the Ising model (for previous studies of the fermion 
operator in this particular case, see [5l[7l[9]). 

4.1 Additional equation 

Let us start by noticing that, in general, eq. ()3.7p is not sufficient to determine the function 
{'ipp{e)). The reason is that, on a covering lattice with N sites, we have to determine the 
complex function {tpp{e)) on 2N edges, but we only have the single complex equation (j3.7p for 
each of the 0{N) plaquettes. 

However, in the particular case Q = 2 and p = ^ it is possible to see that another complex 
relation holds for every elementary plaquette. This has the form 

(V'i/2(ei)) - (V'i/2(e2)) + (^'1/2(63)) - (-01/2(64)) = , 

where ei, 64 are the four consecutive edges of the plaquette. 

In order to formulate this in a more systematic fashion, let us consider again the eight classes 
of configuration in Fig.O Let us define 
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Yl e-*^'2™p(G) , (4.1) 



W = — 00 Q^(j('^) 



where C^^^ is the set of all graphs in the FK expansion which produce configurations of type 
{ih) {i = 1, ..,4), with winding number w, and P{G) is the weight of a graph G. It follows from 
()2.4p that for every G G Gih there is a corresponding G G Gia such that 



P{Gla)=uP{Gu) , P{G2a) = - P{G2b) , P^G^a) = - P{G3b) , P{G^a) = U P{G4b) 

u u 
(at criticality all factors reduce to y/Q). We therefore have 

(V'p(ei)) = |(n+ 1)51 + + 1) ^2 + <73 + 54} (4.2) 
(^p(e2)) = I (n + 1) e^^i 51 + e'"^ 52 + (| + l) e^"^ 53 + e'^^i" 94 

(^p(e4)) = je-'Pisi + (f + 1) ^"'''^52 + e^^i'^53 + (n + 1) e'^^'t ^4 
which can be written in the matrix notation 

ii;) = M-g. (4.3) 
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It is now easy to check that 

14 for generic -^/Q 
3 forVQ = 2sinpf . (4.4) 

2 for Q = 2 , p = i 

Thus in general, if we choose p correctly, there is one linear relation between the parafermions 
around the 4 edges of each plaquette, while for Q = 2 there are two. 

4.2 Continuum limit of the lattice equations 

In this section we make some remarks on the question of the extent to which our results on 
lattice holomorphicity imply that in the continuum limit the parafermion correlators are truly 
holomorphic functions of their arguments. In general this does not follow without further as- 
sumptions. Recall Morera's theorem[8] which states that if ^fjf{z)dz vanishes for every closed 
contour C in some domain T>, and f{z) is continuous, then f{z) is a complex analytic function in 
T>. Although our lattice result shows that there is a subsequence of approximations to f^, f{z)dz 
which do vanish identically, this does not imply that the limit vanishes without further smooth- 
ness assumptions. Although we expect these to be valid, it would require further work to show 
that they follow from the lattice model. Nevertheless, if true, this implies that in, for example, 
the upper half plane when the curve runs from to infinity, {ipp{z)) const. z~^. 

In the case of the Ising model, Q = 2, the situation is better. In that case there are sufficient 
equations to determine the lattice function (^(e)). Since they are linear and translationally 
invariant, it should be a simple matter to solve then by Fourier analysis in a sufficiently simple 
domain, and then to show that, in the limit when the size of the domain is much larger than the 
lattice spacing, the Fourier components near A; = dominate, and in the limit, are given by the 
solution of the continuum problem. However, this approach is not adequate for more irregular 
domains. 

4.3 OfF-criticality 

As soon as we move away from the critical point, eq. (j3.7|) does not hold anymore, since the 
necessary cancellation mechanism relies on the particular value Uc = y/Q- However, linear 
relations between the values of (ipp) at different edges are still present if ^/Q = 2sinp^, since 
()4.4p holds for any value of u. For instance, if we define v = u — \fQ we have 

(^p(ei)) + z(^p(e2)) - (Vp(e3)) - ik^Ae^)) = X {{^.p{e^)) - (Vp(ei))) , (4.5) 

where 

1 + ie'P^ 

X = V —. = . 

gipTT _i_ 2sinp| - V 

We will now give a precise meaning to this off-critical relation in the case Q = 2. It is well 
known that the Ising model can be described, in the scaling regime near its phase transition, by 
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a field theory of a free Majorana fermion of chiral components X) with the action 

5 = ^ y d^z {xdzX + xdzX + imxx) , (4.6) 

where the parameter m measures the deviation from the critical temperature. In the above 
notation, m is positive for T < Tc and negative for T > Tc- The absolute value \m\ is the mass 
of the elementary excitation in both regimes. The action (j4.6p implies the equations of motion 

Tfl Tfl 

dzX = i^X: dzX = -i^ X ■ (4.7) 

The left hand side of (j4.5p is the discrete version of 2idz{'ipp), therefore we have to reexpress 
the right hand side in terms of {^pp) in order to reproduce a discrete version of (j4.7|) . This is 
possible only in the case Q = 2, because at the value p = ^ the functions gi defined in (14. ip 
coincide with their complex conjugates g^. In more general cases, the right hand side of (j4.5p 
cannot be just written in terms of {ipp), and this is consistent with the fact that no simple 
equation like (|4.7p holds in the corresponding field theory descriptions. Going back to the Ising 
case, we have 

(V'i/2(ei)) +^(^1/2(62)) - (^1/2(63)) - «(V'i/2(e4)) = 

V _ _ _ _ 

= , /TTx/o AT , ^ ((V'i/2(ei)) + (^^1/2(62)) + (^1/2(63)) + (V'i/2(e4))) , 
(1 + V2)(2v2 + V) 

which is the discrete version of 

^^^--^^^/^^= \l + V2)(2V2 + .) '^^"^/^^ 

where the antiholomorphic operator in the right hand side is uniformly distributed on the pla- 
quette. This relation and its complex conjugate are equivalent to (j4.7p if the fields are related 
by 

X = aV'1/2 , X = (3^1/2 , with - = i 



and the parameter m is given by 



e 



J 



m 



1 - V2 



^1 + ^2) (2^2 + 1;) (l + \/2)(\/2-l + e- 



V2{J-Je)+0{J-Jc)^. (4.. 



We can now compare this result with the known value of the correlation length ^ on the square 
lattice, expressed in terms of the modular parameter k = (sinh J)~^ [10] and evaluated in the 
scaling limit: 



-logk ^2V2{J-Jc)+0{J-Jcf for T<Tc 

(4.C 

Uogk ^ V2{Jc- J) + 0{Jc- J)^ for T>Tc 



From the definition of the correlation length in terms of the connected two-point function of 
the order operator a 

{a{x)a{0)) ~ x~'^ e~^^^ as x — > 00 , 
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it follows that 1/^ must be equal to the mass of the lightest excitation which couples to a. 
The Z2 symmetry of the Ising model implies that a only couples to states with odd numbers 
of particles at high temperature and with even number of particles at low temperature. Hence 
()4.8p and (|4.9p are in perfect agreement. 

5 Continuum predictions 

In this Section, we will show that, despite difficulties in proving convergence to the continuum 
limit, the lattice holomorphic objects defined above nevertheless have their counterparts in the 
continuum formulations of these models. In particular, we will review results obtained by CFT 
and Coulomb gas methods, and we will describe a new prediction from SLE. 

5.1 CFT and Coulomb gas results 

From the CFT point of view, the holomorphicity condition p.3p has a clear meaning both for 
Q = 2 and Q = 3, where the values P = \ and P = | coincide with the spins of the parafermionic 
currents predicted by Fateev and Zamolodchikov [1] for CFT with extended Z2 and Z3 symmetry, 
respectively. 

A more general result was obtained in [6] for the Potts model by using Coulomb gas methods. 
Nienhuis and Knops showed that for every Q, there is a series of spinor operators parameterized 
by a rational number p < 1, whose denominator is the order of one of the cyclic permutations 
into which the permutations of {0, 1, Q — 1} can be decomposed. For instance, at Q = 2 we 
only have p = 1/2, while at Q = 3 we have p = 1/2, 1/3 or 2/3. The scaling dimension Xp and 
spin Sp of the spinors are given by 

Xp = I + (P^ - 1) , Sp = p , (5.1) 

^ y 

where y is defined through 

= 2 cos (ti I) . (5.2) 

For every Q, there is a single holomorphic operator among these, characterized by Xp* = Sp*, 
i.e. 

p* = l-y. (5.3) 

This result is in perfect agreement with our requirement (13. 3p for the holomorphicity of the 
parafermionic correlation functions on the lattice. 

5.2 The SLE prediction 

In this Section, we will show that some local observables associated to an SLE curve are char- 
acterized by the scaling dimension of the spinor operators discussed above. To this purpose, 
we will recall and generalize results obtained in [3] for the analysis of another operator, the 
stress-energy tensor. 
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Let us consider a chordal SLE^ process (for < k < 8) on the upper half plane described by 
complex coordinates w,-!]) [TT]. The probability P{wi,W2,wi,'W2) that the curve passes between 
two points wi , 'W2 (or to the left or right of both) satisfies the equation 

^ {dwi + du,-, + + dwif + —dwi + -^d'wi + —dw2 + -^9^2 \ Piwi,W2,wi,W2) = . 

Z Wl Wl W2 W2 J 

(5.4) 

If we parameterize the event by the middle point w oi a straight segment, by its length e and 
by the angle 9 that it makes with the positive imaginary direction: 

6 7(9 ^ iO 

Wl = w - -e , W2 = w + -e , 
we obtain, at leading order in e, the following equation 

'^{du, + du,f + -du,+ '^du,-(\ + ^ ed, + ( \ - ^\ ide + 0{e^) \ P{w,w,e,e) = 
2 w w xw^ J yw-^ J J 

(5.5) 

In the above expansion we have assumed that each of the Fourier modes of the probability, 
defined as 

oo 

Qp{w,w,e)= j dee~'P^ P{w,w,€,e) , (5.6) 

— oo 

vanishes with a power law e^p as e ^ 0, hence = 0{e~^). The parameter p has the physical 
meaning of " spin" and took the value p = 2 in the case of the stress-energy tensor studied in [3] . 
Here we are interested in more general values of p, in particular non-integer values. Therefore, 
the domain of integration for the variable 9 in (j5.6|) is not restricted between and 2tt, since 
the winding of the SLE curve around the segment is associated to a shift 6 + 27rn, where n G Z 
is the winding number. If p is integer this phenomenon is irrelevant, but ii p < 1 winding has 
to be taken into account. In particular, there are d inequivalent winding classes, where d is the 
denominator of p. We can define the probability that the SLE curve passes between the ending 
points of the segment with a winding number n = k (mod d) as 

oo 

Pk{w,iB,e,9) = ^ P{w,w,e,9 + 2TTk + 2iTmd) . 



m=—oo 



In this notation, the Fourier modes can be expressed as 

J 1 2n 



Qp{w,w,e) = ^e-^v-i--^^ f d9 e'^P^ Pk{w,iv,e,9) . (5.7) 
fc=0 { 

By integrating eq. (|5.5p over f d9 e~^^^ we obtain, to leading order in e, 



^((9iu + dii,f +'^_dyj+'^ dii, - + ed^-p ( ^ - \ Qp{w,w,e)+coiciections = 

(5.8) 



W W XW^W^ Xw^w^ 
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Among the several solutions of (jS.Sp , we are now interested in the one which satisfies appropriate 
boundary conditions for the event of passing in between the two initial points. This particular 
solution was obtained in [3] by mapping the problem into the disk geometry; for completeness, 
we report the detailed discussion in AppendixlAl The solution is 

Qp{w,w,€) = Cp e'^p w^p {w - , (5.9) 

with 

K-8 p K-8 p (8-k)2-kV -, 1^ K 2 ^^^ 

a„ = , Pn = — , 7n = , x„ = 1 — p . 5.1U 

^ 2k 2 '^^ 2k 2 ' 8k ^ 8 8^' 

The function (j5.9|) gives the correct leading behaviour for small e of Qp{w,w,€) only when 
the higher order terms in (j5.5p do not mix the Fourier components at leading order. This is 
guaranteed if 

Xp < 2m + Xp-2m V ?n = 1, 2, 3, ... such that Cp_2m 7^ , (5-11) 

where m labels the powers of in the expansion of the differential equation (j5.4p . This auto- 
matically holds for any < p < 1, which is the case of interest in this paper. 
At the particular value 

Kp = ^ , (5.12) 

(|5.9p simplifies to the purely holomorphic function 

Qp{w^w,e) = const x {—] , (5.13) 

where spin and scaling dimension are equal. 

In order to compare the SLE prediction with the other results described in this paper, let us 
recall that the Q-states Potts model in the FK representation is conjectured [12] to be described, 
in the continuum limit, by SLE^ where 

g = -2cos ( TT- ) , 4<K<8. (5.14) 



It is then easy to check that (jS.lOp perfectly agrees with the Coulomb gas results ()5.ip for every 
< p < 1, and in particular at the holomorphic value (j5.12p . 

The holomorphic solution (j5.13p has a natural CFT interpretation as the expectation value 
of a purely holomorphic operator of spin p, given the appropriate boundary conditions which 
generate the domain wall. By recalling the relation between k and the central charge c of CFT 

m 

2/c 

one can see that (|5.12p implies 

/i3 1 for < K < 4 , ^ 

P = <; ' , (5.16) 

/ii 3 for 4 < K < 8 



13 



where hij are the conformal weights of primary operators in the Kac table of CFT. For Ising 
and 3-state Potts models, this is again in agreement with the other results presented in this 
paper. Relation (j5.16p . however, is more general and deserves further study in connection with 
other statistical models. 

6 Relation between SLE and lattice results 

In the last Section, we have checked that the SLE prediction agrees with other continuum limit 
results. Here, we will discuss the explicit relation between the purely holomorphic function 
(j5.13p obtained from SLE and the lattice spinor ipp defined in (j3.2|) . 

To this purpose, we have to look at the one-point function {ipp{e)), which is an observable of 
the domain wall that is conjectured to converge to SLE in the continuum limit. We have shown 
in Sect. 13. 21 that (V'p(e)) is purely holomorphic when ^/Q = 2sinp|, i.e. when k = as in 
(|5.12p . Moreover, we know that (^p(e)) satisfies the same boundary conditions obeyed by 
in the continuum limit, hence we expect 

Qp{w,w,e) = const x {ip.p{w)Y'''^^ , (6.1) 

where {'4^p{w)Y°'^^ is the function to which {'ipp{e)) is supposed to converge in the continuum 
limili. 

The advantage of SLE, with respect to other continuum limit approaches, is to have a direct 
and explicit meaning in terms of the lattice formulation of the problem. This allows us to 
analyse the correspondence between Qp{w,'w, e) and {ipp{e)) on the lattice itself, by defining a 
discretized version Qjf* of Qp where the SLE probability distributions is replaced by the lattice 
weights in the FK representation of the Potts model. In principle, one should check the relation 
for segments connecting two arbitrary points wi,W2 on the lattice. As a matter of fact, the 
computation becomes very complicated as soon as several plaquettes are involved, and we have 
been able to carry it out explicitly only for the case of a single plaquette. This corresponds to 
the rotating segment shown in Fig. 21 



(1) (2) (3) (4) 




0=7c/4 e=37c/4 0=57c/4 Q=7n/4 



Figure 4: Rotating segment within a single plaquette. 
■^For p 7^ 1/2, we need a continuity assumption to justify this statement, as discussed in Sect.|ll 
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Let us first obtain the lattice version of Qp. Tlie range of integration in (j5.7p can be further 
reduced to G [0, tt] by noticing that, in every winding sector, the curve can pass between the 
end points of the segment by leaving point wi to its left and W2 to its right, or viceversa, and 
the two events are related by the transformation 9 ^ 9 + n: 



Qp{w,W,€) 



d9 e 



-ip t 



[Pk,iLM^^w,e,9) + e 'P^Pk,iR,2L{^^w,e,9 + 7:)] . 



This implies that we only have to analyse cases (1) and (2) in Fig.[31 The discretisation of the 
integral in 9 gives 



Qp 



lat 



^ k 



^kAL,2R ^ ^ ^k,lR,2L 



,(1) 



-ip Ivr 



^kAL,2R ^k,m,2L 



-ip TT ip 2-K 



+ e 



-ip\-K 



} 



where the functions gi are defined in (|4.1|) . and their relation to the probabilities of passing 
between the end points of the segment can be inferred from Fig.[2j The factors e*^^'^ in the last 
line are due to the fact that the domain wall in configurations (la) and (3a) of Fig. [2] passes 
between the end points of the segment with winding —1. 

We expect from (j6.ip that (|6.2|) is proportional to {'4)p{w)) evaluated at the middle point w 
of the segment, i.e. at the center of the plaquette. Assuming the spinor operator to be uniformly 
distributed on the plaquette, we have 

1 1 ^ 

(^pM) = 4 ((^p(ei)) + (^^^(62)) + (V'p(e3)) + (^p(e4))) = ^Xl^^^^' 

i=l 



where the coefficients Si can be easily obtained from (j4.2p at u 
that 

iiat _ \/Q secpf 



Q. It is now easy to check 



(6.3) 



2 [^ + 2cospf] 
which confirms (j6.ip . 

It is worth noting that ()6.3p is valid for any value of Q, not only for the purely holomorphic 
case \/Q = 2sinp^. This agrees with the general SLE prediction (j5.9p that, for every Q, the 
several spinor operators obtained by Coulomb gas methods in [6] are local observables of the 
domain wall. 

Let us conclude this Section by analysing the slightly different situation in which Qp is 
defined as the Fourier component of the probability that the domain wall intersects the segment 
without necessarily passing trough its end points. Additional configurations have to be taken 
into account (see Fig.[2|), and (|6.2p is modified as 



lat 



Qr + 1 



-tp - 



[9i+- 



-ipn 



which still implies a proportionality with {'il)p{w)): 



54 +6 



" [92 + 



lat 



2tt 



seep J 



2 [VQ + 2cospf] 



(ippiw)) . 



(6.4) 



(6.5) 
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This is consistent with the SLE result, because the boundary conditions that select solution 
(|5.9|) are the same if we consider the events of intersecting the segment or passing between its 
end points (in the disk geometry defined in Appendix Rl both events require the curve to pass 
by the center of the disk) . Therefore, if we assume that the distortion of the segment under the 
Loewner conformal map only introduces terms of higher order in e in (jS.Sp and (j5.8p . Qp and 
Qp are expected to only differ by a constant in the continuum limit, and this is nicely confirmed 
by the lattice example analysed here. 

7 Conclusions 

In this paper, we have carried out a detailed study of parafermionic operators in the Q-state 
Potts model, both on the lattice and in the continuum limit. 

Our discrete setting was a planar square lattice. By geometrically defining the parafermions 
in terms of the FK representation of the model, we obtained an explicit proof of lattice holomor- 
phicity for certain values of the spin. We believe that this is first time that these parafermions 
have been identified within the FK representation of the Ising and Potts models, and that their 
holomorphicity has been directly demonstrated at this level. It would be interesting to extend 
our analysis to different lattices and geometries. 

After reviewing known results on the scaling dimensions of parafermions from CFT and 
Coulomb gas methods in the continuum limit, we have shown that the same scaling dimensions 
are associated to local observables of SLE^ when (j5.14p holds, i.e. at the value of k which is 
conjectured to describe the Q-state Potts model in the FK representation. Moreover, we have 
explicitly discussed the relation between the SLE observable and the lattice parafermions, both 
in the discrete and continuum settings. 

Natural extensions of this work are the full analysis of the situation when the parafermionic 
operators are not purely holomorphic, and the study of the SLE prediction in the context of 
other statistical systems. In particular, it will be interesting to exploit the information provided 
by SLE for models associated with non-unitary or non-minimal CFT. 
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A SLE probabilities in the disk geometry 



The ansatz 

Qp{w,w,e) = Cpe""^ w""^ wf^^ {w -wf^ (A.l) 
solves eq. (|5.8p for two different choices of the parameters: 

/c — 4 At — 4 (k — 4)^ (k — 4)^ 

and 

"^ = ^"2' ^- = ^ + 2' > = ' "- = '-8 + 8^- ^^-'^ 

In order to select the correct set of parameters, it is convenient to map our problem onto the 
unit disk D, through the transformation z' = for z £ M and z' G D. This transformation 
maps the point w to the center of the disk, the length e to e/lw — w\, and it shifts the angle 6 by 
an angle of 7r/2. Also, the point is mapped to w/w on the boundary of the disk, and the point 
00 to 1. We are then describing an SLE curve on the unit disk started at w/w and required to 
end at 1. Fixing the power of e/lw — iDl to be some number Xp (the "scaling dimension"), we 
are left, after integration over 9 as in (j5.7p . with a second order ordinary differential equation in 
the angle a = a.rg{w/id) G [0, 2-k]. This equation is the eigenvalue equation for an eigenfunction 
of the two-particle Calogero-Sutherland Hamiltonian with eigenvalue (energy) 2xp/K and with 
total momentum p |14] . For generic k, the Calogero-Sutherland Hamiltonian admits only two 
types of series expansions Ca'^[[a^]] (with C 7^ 0) as a — > 0"^ for its eigenfunctions: one with a 
leading power a; = 8/k — 1, the other with a leading power u> = 0. It admits the same two types 
of series expansions C'{2tt — a)^' [[(27r — a)^]] (with C" 7^ 0) as a — > 27r~. Allowing only one type 
of series expansion at and only one at 27r (the possibilities give the Calogero-Sutherland system 
in the fermionic sector lj = uj' = 8/k — 1, bosonic sector lo = u>' = or mixed sector, uj 7^ lv'), the 
Calogero-Sutherland Hamiltonian has a discrete set of eigenfunctions, with eigenvalues bounded 
from below (since it is a self-adjoint operator on the space of functions with these asymptotic 
conditions). The lowest eigenvalue is obtained for the eigenfunction (the ground state) with the 
least number of nodes (zeros of the eigenfunction). If the leading powers to and uj' are chosen 
equal to each other, then the ground state (in the sector with total momentum n) is described by 
the solutions (|AT]) with (|Ar2|) (for w = 0) or (for w = 8/k - 1), which, in the coordinates 

of the disk, take the form 

^ n -I N ~ f € \'^'' i °-P-0P ^ ( . a\^v+^v ,. ^, 

Qp(\w-w\,a,e) = Cp\- 2 sm - . (A.4) 

\\w — w\) \ 2j 

If we consider the probability that the curve passes in between the two points in the original 
formulation on the half plane, then the curve is required here to pass by the center of the disk. 
Hence, the probability vanishes when the starting point of the SLE curve is brought toward its 
ending point on the disk, from any direction; this fixes the power to be 8/k; — 1 (for k < 8) 
at both values a = 0,27r and therefore selects the solution in the fermionic sector (|A.3p . Note 



17 



that since the probabihty could be given by an excited state in the fermionic sector (which 
corresponds to a higher value in place of the exponent Xp), we do not have the condition that 
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